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ABSTRACT 



An investigation is made of the wave forecasting method of Pierson, 
Heumann and James of New York University, published as H, 0« Pub« Ho. 603. 
The oooan "swell” spectral energy, as predicted by, the application of this 
forecast method is compared to the observed spectral energy present as 
estimated from ocean wave recording devices installed at the svirfacb and 
at a depth of 80 feet, in the coastal waters off Davenport, California. 

The comparison of the obseinred and forecast spectral energy for nine 
verification times taken at t\relve hour intervals from 28 Janii^y 1954 to 
1 February 1954 suggested the following conclusions; 

1. For the conditions of the investigation, it was the opinion of the 

authors that the "swell" forecast was operationally useful. ^ 

2. The data of this investigation suggested the possibility that the 
forecast energy speotrvaa is too high in the low and middle frequencies and 
too low in the high frequency bands for high wind speeds* 

3. The data of this investigation indicated that at certain verifica- 
tion times, energy reductions due to cross-sea interference were important. 

The authors wish to express their appreciation for the advice and 
encoui'agement of Professor J. B. Wickham under. vdiose direction this study 
was curried out. 

Appreciation 'is also expressed to Professor R. L. Wiegel, of the 
Institute of Engineering Research, Berkeley, California, who made 
available the wave records upon which this study is based* 
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CHAPTER 1 



DITRODUCTION 

The value of reliable techniques for forecasting ocean vmve proper- 
ties for seagoing and coastal operations, particularly military operations, 
is obvious* 

This problem has engaged the interest and attention of some of the 
foremost scientists of various fields over the years* Somo of tho contri- 
butions have been notable. 

The theories of classical hydrodynamics were applied by Lamb Cl] 
in 1879, which yielded important properties of certain relatively simple 
sine -wave forms* 

In 1933, Comish[^23, published the results of a large number of 
carefully made visual observations of wind-vraves, both in deep v/ator and 
near the shore, which wore the first of that kind and vdiich undoubtedly 
stimulated further work of this nature. 

The development of continuous wave recording devices expanded the 
vmve data available and made possible the many statistical studies which 
have been carried out, 

Sverdrup and Hunk [4j published a wave forecasting method in 1947, 
for predicting tho generation, decay, and transformation of the 
"significant” wind-waves* 

Bret Schneider, in 1951, published a modification forecast [sj of 
the Sverdrup and Hunk technique based on more observational data. 

In 1953, Pierson, Neumann and James of New York University pub- 
lished a new theory and corresponding forecast technique * This 
theory and method are in many respects unique and depend largely on some 
assumptions vdiich have had limited testing. 
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Sinoe this last theory and toohnique has been adopted for use by the 
United States Naval Service, a study and evaluation of the nothod is the 
purpose of this paper* A study of this nature presents many difficulties 
and requires data and facilities available at very fevv institutions* These 
facts explain the paucity of published results of such studies* To the 
knowledge of the authors, no such studies and evaluations havo been made 
previously of the Pierson, Neumann and James wave forecast method, hence- 
forth referred to as the N. Y, U. method, for Pacific Ocean swell, although 
the waves forecast in a generating area havo been studied [isj. 

The study divides itself naturally into four distinct phases: 

1* A wfivo forecast for a given period of time at a specified location, 
using the N. Y, U. method, 

2* The analysis of wave records from a wave recording device for the 
same location and same time to obtain the actual vrave conditions* 

3* The comparison of the forecast and observed wave properties. 

4* The results or conclusions to be dravm from the oemparison* 

The data and facilities necessary to carry out this study are listed 
below vdth their sources: 

1. A forecast location at vdiich suitable v/ave records have been 
kept for a convenient period of time. The site chosen was Davenport, 
California, Tirhero the University of California, Berkoloy, had a series of 
recording devices LO installed from the period of Nov* 1952 to Dec* 1954. 
\loxc records from Jan. 1954 to Feb* 1954 were made available to the 
authors for use during this study* 

2* A rapid computing device for the analysis of the wave records 
to obtain the distribution of wave energy as, a function of frequency [V). 

The U* S* Naval Postgraduate School is equipped with a CRC 102-A medium- 
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speed c®ii6ral-purpos© electronic computer which was available for this purposso 

3. Synoptic surface weather maps for the period of and preceding 
the forecast time for the purpose of preparing the detailed wave forecast. 

Such maps vrere made available by the Department ef lloohanioal Engineering 
V/ave Research Laboratory, University of California, Berkeley. 

4, A submarine contour chart for the area adjacent to the forecast 
location. U, S. Coast and Geodetic- Original Survey Sheet llo. 5266 of 

25 April 1935, was made available by the Hopkins Marine Station, Pacific 
Grove, California. 



3 



CHAPTER 2 



THE THEORY Al'ID THE LIETHOD 



The objective of this chapter is to discuss the portions of the 
H. y. U» Method and theory -which relate directly to this in-vestiga- 
tion. 



EA The Spectrum of Wind- Gene rated Waves 

The primary basis for this forecast method is the assumed distribu- 
tion of ocean vrave energy for a given vdnd speed as a function of -wave 
frequency in a ■wind--wave generation sirea, or the energy spectrnmi. This 
distribution is represented mathematically by: 



li%ere c is a constant of proportionality and the other terms are defined 
on page VI* Figure 1 shows a plot of this function for -wind speeds of 20, 



theoretical considerations and empirical evidence* To a large measure, 
the aooxiracy of this method will be determined by how closely the assumed 
spectrum approximates the true spectrum* Such determination -will be 
possible only as increased test results became available* 

2A1 Energy in the Spectrum 

The integral of Equation (l) for a gi-ron wind -srelocity vrith respect 
to frequency is proportional to the energy of the included spectral com- 
ponents* This measure of spectral energy, bo-bween the frequencies 
and ^ 2 , * is defined as E and is given by 




( 1 ) 



30, and 40 knots, and tabulated values of 




function of wind 



speed and frequency are given in Appendix I* 

The form of this function was assianed by Neumann (jLO^I based on 
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df 



( 2 ) 




and has the dimensions, length sqp^ared. 

Well-tested statistical relationships [ll*] between E and v/ave heights 
haTo been demonstrated. Those relationships are given in Table 1 and are 
used as a basis for all wave height dotorminations, E being tho forecasted 
parameter* 

2AZ Properties of tho Spectrum 

The wind-generated sea for a given mnd speed is said to be fully 
developed vdien all possible components in tho spectirum botvraen f = 0 

and f s oo are present with their maximum amounts of spectral energy. 

This state is reached only after the wind has blown for a sufficient 
length of time C’diuration” ) and over a sufficiently long distance ("fetch"). 
If the wind duration or fetch is insufficient for the fully developed 
state, a lower energy level will result. The spectral components form in 
the high frequency end of the spectrum initially and the lower frequency 
components are added progressively and grow as the wind dxxration increases 
until full development occurs, providing the fetch is not limiting and the 
\Tind continues. One important feature of this theory shows tho partially 
developed spectra for limiting conditions of fetch and d\iration for 
various wind speeds. 

EA.3 The Co-cumulative Spectrum 

Tho energy spectrum is converted to a more convenient form by 
integrating it ovor tho entire frequency range. This procedure is illus- 
trated in Figure 2. It vrill be noted that the ordinate of the cixmulative 
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integral can be read directly in values of E, which is proportional to all 
of tho ivavo energy in spectral oonipononto equal to or greater than tho 
abscissa, f. The functional relationships betirroon energy, v/ind, duration 
and fetch may also be included on tho cumulative spectrvim. 

Figure 3 and Figure 4 illustrate the 1I« Y, U. co-cunulativc spectra 
given as functions of duration and fetch respectively for various v^ind 
speeds* These are the cxirvos normally used in tho II. Y. (J. method. 

2B Tho Wave Forecast in tho Generating Area 

The foregoing are required for forecasting the wave conditions in a 
wind-wave generating area. Kna^ring the vrind speed, tho wind duration, and 
the fetch length, the E-value and the range of wave periods with import- 
ant energy may be obtained diroctly from tho co-cvmulative spectrum. Tho 
sea in the generating area has an irreguleir, short-crested, appearance as 
tho spoctral components randomly amplify and damp each other. The 
spectral components also move in a variety of directions, as the \Tind 

I 

action is not uniform, but of varying strength and direction. 

20 The ’’Swell” Forecast 

VIhen wind-gonerated waves leave the generating area, a transforma- 
tion takes place in the appearance of the sea surface as tho wind no 
longer acts on tho wave form. As the various spoctral componont vmve 
trains leave tho vrind area, the apparent vmves become regular and long 
crested, and arc termed "swell.” lYhon "swell" arrives at points distant 
from the generating area, tho v/avos have only a fraction of the energy 
originally present in the fetch. 

201 Assumed Processes Involved in Energy Reductions 

Pierson, lleunann and Janos accovint for energy reductions at points 
distant from tho fetch by tiro processes, namely 
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1. Dispersion 

2* Angular Spreading 
202 Dispersion 

When the spectral oomponents leave a generating area, the associated 
energy front of each travels at different group^ velocity. The group 
velocity for any spectral frequency is given by 

G. V. a 1.515 T ■ (3) 

Dispersion, then, is defined to be "the spreading-out effect caused by the 
different group velocities of the spectral frequencies in the original 
disturbance at the source." (6j Therefore, to determine for a given time 
the total spectral energy at a specified distance from the fetch, based on 
the effects of dispersion, it is necessary to compute the range of fre- 
quencies vdiich could be present at that time. This may be accomplished by 
considering the distance the energy must travel, the group velocity, and 
the time the particular spectral component was formed as a function of 
wind speed, duratien, and fetch. 

The IJ. Y. U. method includes some basic "filters" "which are essen- 
tially pre-arranged formulations of these "vnriables, for certain recurring 
situations. The use of those "filters" permits the calculation of a T^, 
upper period, and Tj^, lawer period. T^ and T^ define the upper and lower 
limits of the spectral energy present at the forecast time and location, 
for a given wind speed, 

2G3 Angular Spreading 

The considerations of paragraph 2C2 do not taJeo into account the 
effect of variability in the direction of wave motion and the fact that 
the "waves are short-crested. The component wave trains vaiy not only in 
frequency but also in direction within tho area of generation. Thus, as 
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tho spectral oompononts leave the generating area, they spread out 
angularly as -well as ’’disperse." 

The assumed distribution of relative wave energy as a function of 
direction is illustrated in the lower curve of Figure 5. Since the area 
under this ctirve is unity, vdien integrated cumulatively over direction it 
gives that fraction of tho energy to be found botvreen any tiTo directions 
(measured from that of tho mean wind). The ordinate of the upper curve in 
Figure 5 is then expressable in percent. Dispersion and angular spreading 
computations are illustrated in chapter 3.”" 

Figvire 6 illustrates the measurement of the angles and O3, vdiich 
are the parameters used in conjimction with the upper curve of Figure 5 
for evaluating tho effect of angular spreading. Sie convention adopted 
for measuring and ©3 is that clookvrise angles are positive and counter- 
clockvfise angles are negative* 

The difference between the percentages obtained by entering the upper 
curve of Figure 5 with and ©3 is tho percentage within that range of 
angles of the spectral energy computed from the dispersion envelope and 
the energy spectrum. 

The bell-shaped distribution of energy with direction. Figure 5, was 
inferred from limited evidence [12-] and to tho knowledge of the authors, 
its exact form has had limited observational verification. 

204 Other Sources of Possible Energy Loss 

Pierson, ITeumann and James have stated that in many cases, 
dispersion and angular spreading are the only mechanisms necessary to 
account for the energy decrease observed in vmve spectra at locations 
remote from tho generating area. They held that there may be important 
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onorgy lost due to the effect of cross-sea interforonoo,^ but have 
expressed the opinion that losses duo to viscosity are generally negli- 
gible. The effects of cross and opposing Virinds on wind-generated swell 
arc not knovm, but undoubtedly would modify the spoctrvm of the swell to 
an extent dependent on the wind speed and the spectral composition of the 
cross-sea. It is believed that at least a qualitative evaluation of these 
effects v/ill bo possible as the volume of test results boo cane greater. 

The direct measurement of attenuation presents serious difficulties. 
Hovrevor, a careful analysis of the actual spectral energy at a wave 
recorder as compared to the forecast spectral energy may indicate the ' 
presence of processes which cause energy losses. For instance, theory 
indicates that viscosity effects, if they eire important at all, reduce the 
energy primarily in the high frequency end of the spectrum If high 

frequency energy freon a distant fetch is forecast to arrive at a wave 
recorder at a given time, and fails to arrive, it might be inferred that 
viscosity effects wore responsible for the energy loss. This possibility 
is discussed further in Chapter 6. 

2C5 Goals of the Investigation 

The goal of this paper is to make available the results of the 
comparison of the spectral onorgy forecast by the N. Y. U. method and the 
spectral energy as obtained from a wave-recording device. 

The assumed effect of cross-sea interference is that the swell with 
spectral components of the same frequency as the cross-sca wave trains 
will form unstable interference patterns with the latter. This vri.ll re- 
sult in attenuation of the s;vell having spectral energy in those frequencies. 
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It was intended that these results would permit conclusions parti- 
cularly in the follcairing areas s 

1* The assumed form of the energy spectrum. 

2. The assumed form of the wave-energy directional spectrum. 

3. The effects of viscosity and oross-seas on the spectrum of the 
svrell waves. 

4. A general verification of the system which implicitly includes 
the iteios above. 



/ 
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CHAPTER 3 



THE FORECAST OF THE EWERGY SPECTRUM 

This chapter contains a brief summary of the procedures followed by 
the authors in propeiring a forecast of the energy spectrum for a given 
observation time* 

When a breeze comes up, the sea sxirface instantaneously becomos 
covered with tiny ripples iThich form more or less regular arcs of long 
radius* They increase rapidly in height as energy is transferred to the 
waves by the pushing and dragging forces of the wind* Soon a state of 
irregular composite wave motion, called the sea, exists* The fully 
developed state over a long fetch with a sufficiently long duration 
finally will be attained when the wave- generating traction’s have increased 
tho total wave energy to such a point that dissipation balances the work 
done by pushing and dragging forces on the sea sxnrfaoe* The total amount 
of wave energy accumulated in the composite wave motion is now distributed 
over a wide range of frequencies* The curves and assmptions of the 
H* Y* U. method were used in forecasting the energy spectra in the fetch 
area* H* 0* 603 contains several filters that can be applied to account 
for the reduction in the energy spectrum as its components move from the 
fetch area to an observation point* The authors, after considerable 
invosti-gation and study, concluded that the idealized filters in H, 0. 

603 vTOuld not be convenient for use in this study. To achieve the de- 
sired results, a new method of solving for dispersion and anguleir spread- 
ing was devised* The new filter used when solving for the effects of 
dispersion allowed for wind variations in the fetch and resulted in a 
continuous dispersion envelope of periods present at the forecast point. 
The angular spreading factor allowed for any orientation of the fetch* 
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3A Selection of Data 



Appendix II describes the elaborate instrumentation used by the 
UniTorsity of California to record continuously the water level at 
Davenport, Califomia. An examination of the wave records available 
shOTfod that the average wave heights wore higher in the vrinter at Daven- 
port due to well-developed storms in the North Pacific area* Therefore, 
the months of January, FebnAary, and March 1954 wore selected for closer 
study* 

Weather maps at six-ho-ur intervals were then studied for these three 
months* All possible wind circulations that resulted in fetches acting 
with components to?mrd Davenport were recorded* An envelope showing the 
periods (1^ and Tj^) and their estimated times of arrival (ETA) v;as com- 
puted for each fetch* These data (T-g, T^, ETA) were tabled and plotted.^ 
After detailed examination, the period between 0630Z, 23 January 1954 and 
0430Z, 1 Febmary 1954 was selected for thorough investigation and study* 
This time interval was the most favorable for the follcn'iing reasons: 
1* A stationary nearly uniform fetch prevailed* It is discussed in 
Appendix III and illustrated in Figures 7 and 8* 

2* Representative wave records were available in "fast time"'^ with . 
a time interval of twelve hoAirs, as explained in Appendix II* 

3* Fetches A and contained weather station ship "P*" "P" is 

« 

loannod by experienced aerographers * Therefore, its weather reports vrere 
assumed accurate* "P" also reports sea-water temperature (Tg) ivhioh 

^Figure 9 A is an example of this from 0630Z, 23 January 1954 to 0430Z, 

1 Fobrusiry 1954* 

2"Fast Time" is used to indicate the portions of the record I’diere sea 
level was recorded on the scale of 6 inches per 1 minute* 

^ See Appendix III. 
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could be used in computing the stability factor,^ which is used in 
estimating the wind speed. 

4. Only fetches A, B, and C, as defined in Appendix III and illus- 
trated in Figtn*es 7 and 8, were believed to be producing significant 
energy at Davenport. 

5. Local winds in the vicinity of Davenport vrere blovring from the 
land toward the sea and were li^t. 

6. All instnanents were cleaned and calibrated on 24 January 1954 
and assumed to bo in good vrorking order. 

Considerations 1, 2. and 4 established wave forecast times at 
Davenport which are as follows : . 

ti 0430Z 28 Jan 1954 

tg l^OZ 28 Jon 1954 

tj 0430Z 29 Jan 1954 

t^ 1630Z 29 Jan 1954 

tg 0430Z 30 Jan 1954 '' 

tg 1630Z 30 Jan 1954 

ty 0430Z 31 Jan 1964 

tg 1630Z 31 Jan 1954 

tg 0430Z 1 Feb 1954 

Hereafter, the forecast time virill be referred to as t^. 

3B Compilation of Data 

Table 2 records the data measured from the weather charts illustrated 
in Figure 7. 

^See Table 6. 
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Tablo 3 records the data oomputed frcsa parameters listed in Table 2 



by the following prooeduros {symbols aro all defined in Tables 2 and 3): 

!• Stability Factor (Vgg^g) depends on Tg minus Tg^ and is determined 
from Table 5. 

2* Isobar cvirvature correction (Iqq) depends on isobar curvature (l^) 
and stability factor and is determined from Table 6« 

3* Geostrophio vdjid reduced tc sea surface (V ) 

''^gs r (Vg) (Ice) (stability factor) (4) 

4» Svirface wind speed (V) obtained by comparing Vg-jj, the reported 
surface wind, with Vgg, 

5* Limiting fetch (Fj^)» upper period 0^^)» siid E were read directly 
from Figure 4 using V* 

6. Equivalent duration (Tg) was established by comparing generation 
time of last V to present V, for an equal E-value* 

7* Travel time (Tg-jj) 

R /Ftp - Ft 

( 5 ) 

1.515 Tu 

HVhere T-y leaves the v/ind area from near the windward side of the fetch*^ 

8* Loiver period (Tj^) 



Tl = 



( 6 ) 



1.515 Tg^ 

This Tj^ will have the some ETA as T^ because it left from the leeward side 
of the fetch. 

9. Duration time (T^) 

T^ = Tg/AT (7) 

IThere /b-T = 6 hours, or map interval. 



see Figure 7 
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10. Estimated time of arrival (ETA) 

ETA - map time -f Tob (8) 

3G Development of the Filter 

The values of must bo reduced by filters as the spectral compononts 
move from the generating area to the observation point* Chapter 2 explained 
the theory of filters involved in the N. Y« U, method. Two processes aat 
to limit the energy at the obseirvation points dispersion and angular spread- 
ing* 

3D Dispersion 

Dispersion results from the presence of different group velocities 
among the spectral components in the fetch. The writers, to accurately 
account for the effects of dispersion, plotted T^ and Tj^ as a function of 
time.^ Tu left from the w±ndvrard side of the fetch with 
distance to travel* Tl moved out of the leeward side trith distance Rq to 
travel* Only periods between T^ and Tj^ were theoretically present at 
Davenport because of the disp>ersion influence* 

Fetches A, B, and C were fully^developed wave generating areas; 
consequently, upper period waves (Ty) wore also leaving from the leeiTard 
side with travel distance Rq* Their ETA was computed by 



^ob “ 



R. 



( 9 ) 



1*515 Ty 

and Equation (8), Due to a shorter travel distance this Ty arrived ahead 
of a simileir Ty which left from fetch vondward* 



^Soe Table 3 
^See Table 9A 



16 



To include the first arriving -ntives froni the generating area, the 
■waves created as the ■wind continued and the ivaves leaving the fetch area 
■when the wind ceased. Figure 9A was modified into Figure 9B, a dispersion 
en^velope, A dispersion envelope is a plot of Ty and Tj^ present from the 
same fetch as a function of the same ETA* The envelope was plotted frcoi 
values obtained through the use of Figure 4 and Equations (5), (6), (8), 
and (9)* The results are shovm. in Figure 9B ■which illustrates ■the various 
limiting envelopes of periods due to dispersion affecting Davenport before 
and during the forecast periods from fetches A, B, and C* 

The resulting reduction in E due to dispersion is illustrated in 
Figxire 10 for t]^ and Figure 11 for t^^, 

3E Angular Spreading 

Angular spreading results from the -variability in vrave direction in 
and after leaving the fetch areas The angular spreading factor is related 
to the orientation of the observation point with respect to the generating 
area and is therefore dependent upon the parameters ©gy, ®4 Lp ^^, 7 * 

■v'shioh were measiured from the synoptic -weather chart. For times near 
the^ end of a fetch life, important wave energy may come from the rear of 
the fetch (fetch windward)* At these times the angular spreading factor 
decreases in value, consequently, loss energy -will be forecast to arrive. 
The factors for the fotch leeward and fetch -windward were obtained 
from an angular spreading factor graph* The difference be-fevroon these 
-tvTO feiotors was (from leevrard to windward) as great as 20,o* The authors 
considered the possible error introduced by using an angular spreading 

^See Figure 7 

^See Figure 5 ' 
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factor averaged over the entire fetch as unacceptable. To accvirately 
determine the angular spreading value as a function of distance from the 
fetch front. Equation (lO) vms usod. 

• A D a 1.515 (Tob) T - Ro (lO) 

^ D is the distance in nautical miles measured from the leevrai'd side of 
the fetch into the generating area.^ 

Tob travel time. 

T is any period between Ty and Tl for t^^ as determined by the dispersion 
envelope. 

Rq is the distance from observation point to the center of fetch leevrard. 
In Equation (lO), T was replaced by l/frequency, a was computed 
for each 0.005 ohango in frequency for each t^ involving fetches B, and 
0, The distance, was thon calculated and plotted, and a corresponding 
©2 and measxired. These angles, applied to the angular spreading graph, 
yielded the angular spreading factor. This procedure was usod Tdien it was 
evident that the waves were not leaving from the front of the fetch. 

The resulting reduction in E due to angular spreading is shovm in 
Figure 10 for t^ and Figure 11 for t 4 » 

It is considered signifioant that the use of Equation (lO) allcrvred 
for variations in angular spreading due to variations in (RqT^^SlD). 

Figure 12 is a graph of angular spreading factor with respect to AD, for 
fetch A. This factor did not decrease at a constant rate in going from 
fetch leeward to fetch windvrard. Graphs, not included hero, yielded 
similar results for fetches B and C. These variations in rate of change 
of the angular spreading factor depend upon the orientation of the various 
fetches with respect to Davenport. 



^Soo Figaro 7 
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3F Shallow Water Effects 



Shallow water effects of refraction, shoaling, percolation, dif- 
fraction, and friction influenoe E as the deep water wave energy approaches 
the Davenport area* 

Refraction was determined using the method presented in H* 0* 234, 
and a wave refraction diagram constructed using U. S. Coast and Geodetic 
Original Survey Sheet Ko* 5266 of 25 April 1935* The chart shows the 
bottom oontotirs seaward of Davenport to be noeirly straight and parallel. 

The range of deep water wave directions was very small, therefore refraction 
vras treated as if all deep wnter waves eirrived from one angle* Shoaling 
effects were also obtained from H. 0* 234, Plate I* 

Percolation and friction vrare assumed negligible duo to the narrow 
v.'idth of the oontinental shelf. Diffraction was assumed to be negligible* 
All instruments were located in open water vroll away from rook formations 
or man-made stnactures. 

The shallow-water effects of refraction and shoaling wore combined 
into one correction factor. Their influence on E is shovm in Figure 10 
for t]^ and Figure 11 for t 4 * Figures 10 and 11 show the resultant ocean 
wave spectral energy forecast to arrive at the recording instruments for 
t]_ and t^* The area under the oxirve can be integrated by a planimoter 
to evaluate the total energy present for any forecast time* The results 
are shovna in chapter 5, and are plotted with respect to time in 
Figure 14* 
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CHAPTER 4 



SPBCTRU1.I ANALYSIS 

The objective of this chapter is to discuss tho general procedure 
and theory employed in obtaining the energy spectrum from the record of a 
vravo-rooording device* ^ 

In this study, wave records vrero used from both a surface wave record- 
ing device and a pressure recording device.^ The distinction between the 
tv;o is that the former records surface waves and tho latter records wave 
pressures at tho depth in vdiich it is installed* 

4A Theory 

Many of the features of sinrface, and below surface waves have been 
explained by Pierson and Marks and Heumann [l3^ based on the assumption 
that a wave record as a function of time is the result of a Gaussian 
process. The theory vms developed by Tvikey and Hamming among others* 

]h tho publication cited, the latter explain a method vhioh may be used 
for estimating the power spectrum of a steady-state sea surface using the 
Foxa'ior transform of the autocorrelation function* A text vdiioh treats 
correlation Ametions and power spectra in detail is Information TheoryD?3 
by Goldman* 

Tirame and Stinson {l4j present curves of the pressme response 
factor as a function of wave period, T, and depth, H* The pressure 
response .factor is a measure of the attenuation of wave energy vri.th depth* 
4A1 Hvoaerical Estimation of Non-normalized Autocorrelation Function 

and the Pressure Power Spoctrtaa 

The results of Tukey and Hanming are used to derive the power 
^Soe Appendix III 
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spectrum from a pressure record placed at the sea floor in a depth of 
water, II. The possible errors -nhich enter into the numerical analysis 
have been considered. The corrections v/hich must be applied to the raw- 
results and the degree of reliability have boen determined. 

The formulas for niaaerical estimation of the power spectra as 
presented by Tvikey and HammingC?] and as used by Pierson and Larks are 
given by: 



tg - ti, = t3 - t2 = tij - tjj_i s ^t 

P(t^) P(t^.p) (P=0.1.2, m) 



H-P 
^ 



S “ N-P 



n=l 



(11) 



( 12 ) 



m-1 



L = — (Qo/2 > Qp Cos CoslTh) (13) 

^ m ^ m 



P=1 



(h=0,l,— m) 

= .23I^_i/.54 ^/.23 (h*-0,l, m) 



A tm 



|^ph(*^hO 



VhAt 



/m 



(U) 

(15) 

(16) 



T/here/{ r 2 TT f 

Equation (11 ) states that ll-l equally spaced time intervals are 
marked off on the record by points spaced At units apart. The values 
of pressure at these points in terns of departures from the moan are then 
tabulated as P(ti), P(t2) P(tn), — — P(tjj), II numbers resulting. 

The non-normalizod lag coefficients are then foimd by the use of 
Equation (12). ranges from Qq to C^, or m 7^ 1 total mmibors. Qq is 
the Eph meix, or the total spectral energy unoorrected for the scale factor 
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of the record, at the depth II. 

The sumraation of Equation (12) vrae accomplished by the use of a 
CRC-102-A computer.^ 

The ”ravr” estimates of pov'rer in the/^ band from/C< = 

tm 

^ ^ 77_i£jr^ found frcan Equation (13). A plot of the "raur” estimates 

^tm 

with values of h result in a quito jagged graph. The jaggedness is intro- 
duced by the measurement of P at discrete points and by the finite length 
of the record. A smoothing device is given by Equation (14). 

The final number, is an estimate of the moan value of the power 
system from TT ih - ») 

to T r(h V 

At the point given by Equation (15) 

/5k tm ^ tm 

the best estimate of the value of pph is given by Equation (16). 

The result is a smoothing of the "jagged" curve of Equation (13), and this 
continuous c\xrve is the best estimate of the pressure power spectiom. 

For the purposes of this investigation, the computer program for the 
Fourier analysis vms limited to Equation (13). Tho refinements. Equations 
(14 ), (15), (16) were not included in the solution, therefore, the resulting 
spectra are the "raw" or unsmoothed ovirves. This limitation was dictated 
by time considerations as the preparation of such a program is extremely 
longthy. These refinements will be included in the computer program as 
time permits and will be available for similar investigations in the 
futvtre. 

4A2 Reliability of Results 

The final results of paragraph 4A1 are estimates. Estimates can be 
in error, and the important part of this analytical method is that it 
also presents an estimate of the error. Tukey and Hamming ra have shovm 



^hc computer programming of the equations used in this investigation 
vras accomplished by Professors J. B. Wickham and H. H. Llartinez of the 
U. S. Naval Postgraduate School. 
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that tho errors follow distribution v^ith f degrees of freodon given 

by 



f - 



IJ - V m 



(17) 



Givon f, thon the use of Table (7) yields information on tho reliability 
of the estimate for each frequenoy band obtained by the analysis. 

For 20 degrees of freedom, theUj^ obtained from Equation (16 ) may be 
multiplied by 1.8 and by. 63. The true value of U in the band under analysis 
v/ill lio betvreen l.SU^ .63U 9(^3 of the time. 

4B Procedure 

For this investigation, the verification periods were of tho order of 
twenty minutos on the wave record. The At was chosen as two seconds as 
determined by the criterion of Pierson and Marks The total number 

of measurements per vmve record, N, xras of the order of 600. For pro- 
graming convenience, the total nxanbor of lags, m, was taken as 43.^ 

Tho total nvtmber of consocutivo vorifioations, 12 hours apart, was 

2 

nine. It was originally intended that an autocorrelation program and a 
Fourier analysis vrould be carried out for each of those times. However, 
some initial difficulties in tho programming routine restricted the time 
available for the investigation, and time did not permit the fulfillment 
of this objective. 



In future work, the irave record should be analyzed over a longer time 
interval resulting in greater value of li» This will increase the degrees 
of freedom and permit better resolution of the spectral curve. The 
number of lags taken should be increased, 

2 

See pp 13, Chapter 3. 
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CHAPTER 5 



STATELIEirr OF RESULTS 

The purpose of this chapter is to present the results of this 
investigation. 

5A The "Filter" 

It is the opinion of the authors that the "filter" devised for the 
energy spectra forecast in this study represents an accurate, practical 
method of solution. This is particularly true for the case of meirkod 
variation of wind speed \'irith time in a generating area. For example. In 
the case of increasing wind speed with time, energy fronts generated by 
the earlier, lower speed winds will frequently be ovortalcen by energy 
fronts generated by the later, higher speed winds. This additional 
energy, from the higher speed v/ind might be overlooked if a plot of the 
"dispersion envelope" with time were not maintained.^ 

The filters and examples described in the IT. Y. U. method 6 
primarily treat the simplified cases of constant or slov^ly changing wind 
speed. This permits the energy to be obtained directly from the G.C.S. 
curves for one wind speed. 

It must be vinderstood that the "period envelope" must be identified 
precisely v/ith the vfind speed tdiich generated it. It appears that the 
filter herein described^ adequately handles the more complicated 
situations. 

It appears that in forecasting by this method, using the filter 
presented herein, the forecaster will frequently need to vrork directly 
with the spectral curves in addition to the C.C.S, curves. For this 



^See Figure 9B 
^See Chapter 3 
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reason, the tabled solution to Equation (1) is included in Appendix I. 

It is reccmmended that personnel using this method plot the curves of 
the spectra for various wind speeds on an enlarged scale for the purpose 
of representing the spectrum present at the locality and time of interest. 
SB Verification of the Forecast 

Spectrum forecasts were prepared for Davenport, California, for the 
times tj^ throu^ tg as described in Chapter 3. These forecasts were 
prepared for two geographical sites^ having water depths of 50 and 81*5 
feet respectively* A description of these sites and the instirumentation 
of each is presented in Appendix II* 

Data available for verification of these forecasts were the 
following: 

1, Signifiosint wave height measurements for the sea siirface* water 

depth 50 feet, hereafter called site I* These data were computed from 

sea surface wave records by personnel of the University of California, 

Berkeley, and made available to the authors for this study* Significant 

wav© height, H i , is defined as the average of the heigihts of the one- 
Z 

third highest waves* E values may be recovered by the statistical 
relationship. 



A plot of the E values obtained in this manner for the times t]^ through 
tg is illustrated in Figure 14* 

2* A sea-level wave record for site I obtained from a Beach Erosion 
Board, Step-resistor gauge as described in Appendix II* These records 
were analyzed by the procedures described in Chapter 4, and an estimate 

^See Figure 13. 




(18) 
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of tho E value for times ■tl and tg were obtained. 

3. A "pressure” recerd for the site in the water depth of 81.5 feet. 
This site will be referred to as site II. This pressure record was ob- 
tained frcaa a Mark IX pressvire-head instrument described in Appendix II. 
The instnmiont vms installed 1^- feet from the bettom at site II or in a 
irater depth of approximately 80 feet. These records wore analyzed by the 
procedures described in chapter 4 and an estimate of the E-value for 
times t 3 through tp were obtained. 

As described above, the "records" available were obtained from two 
different instruments. For the times and t£p the surface step- 
resistance device was available at site I. For the times tg through tg 
the pressure-head device, installed at a depth of 80 feet, was available 
at site II. 

Fig\are 14 is a plot of E-valuos for the times t]^ through tg for the 
sea surface at site I, as estimated by: 

1. The forecast. 

2. The H 1 computations provided by the University of California, 

S' 

as described above. The confidence limits are plotted at each 
point. 

3. The ivriters' analysis of the water level record of the step- 
resistance gauge as described above. The 90/o confidence limits for 
these points are also included.^ 

Figure 15 is a plot of E-values for tho times tg through tg for 
site II for a depth of 80 feet, as estimated by; 

1. The forecast. 

^See pp 22, Chapter 4. 
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2* The -vjritors' analysis of the "press\ire” record of the instnjmont 
as described above. The 90/t confidence limits for those points are also 
included,^ 

The prooodure for obtaining the forecast spectra and E- values for 
verification times t 3 through tg at site II was the follov/ing, 

A surface spectrum forecast for a given verification vms prepared at site 
II* This spectrum was reduced to a depth of 80 feet using the pressure 
response attenuation factors of Timrae and Stinson Q.43 as described in 
Chapter 4. An example of the results of such a procedure is illustrated 
in Figtires 16 and 17 for verification times tg and tj respectively. 

As stated earlier, the original intention was to obtain the computed 
spectrimi for each verification time, Ilowevor, because of some initial 
difficulties encountered in the computor programs discussed in Chapter 4, 
the time allotted for the use of the computer did not permit the fulfill- 
ment of this objective, 

A discussion of the figures described above, and the conclusions of 
the authors will be presented in Chapter 6* 

^See pp 22, Chapter 4. 



26 



CHAPTER 6 



CONCLUSIONS 

Tho purpose of this chapter is to oompare the obsorvations with the 
forecasts and to suggest possible explanations for any differences noted. 
6A General Comments on the Results 

An examination of Figures 14 and 15, Tdiioh compare forecast and ob- 

r 

served energy, and Figure 22, which compares the forecast and observed 
significant height, indicates certain interesting features. 

1. It is noted that the general agreement is quite good. Tho fore- 
cast energy peaks occurred at the time predicted and the energy was of the 
correct order of magnitude. 

2. The forecast energy was in general equal to or greater than the 
observed. 

3. At the point of greatest difference the forecast significant 
height was about one and one-half feet greater than the observed. This 
would indicate that for these conditions, the forecast would be of 
operational value. 

4. It is noted that the points of closest agreement occur vjhere the 
observed energy is at or near a peak. 

6B Possible Explanations for tho Differences between the Observed and 
Forecast Energy 

For the purpose of sviggesting possible explanations of differences 
between the obsoi^ed and forecast energy, the following items are con- 
sidered to be possible soiirces of error: 

1. Tho theory. The forecast energy depends upon 

a. the energy spectrum and, 

b. the directional spectrum. 
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If the aBsumed form of either of those ouinres is inoorroct, tho energy 
forecast would probably be in error. 

2. Tho application. In applying the method, it is xmavoidablo that 
the selection of data is largely subjective. It is possible therefore that 
sizable errors could bo introduced. Errors in the folloiiVing might be 
particularly critical: 

a. geostrophic vrind analysis 

b. surface wind calculation 

c. measuroiiiont of fetch parometors. 

5. Energy transformation processes. Processes or conditions vdiich 
eire believed to affect ocean swell wave energy and are not occomted for in 
tho forecast are the follcrwing: 

a. oross-soa interference 

b. cross or opposing wind 

c. viscosity effects, 

4. Refraction and shoaling effects. 

6C The Energy Forecast 

Figure 14 compares the forecast value of E vrith the observed at site I 
at the sea surface for tho verification times ^1 through tg. It is noted 
that the forecast E value lies above the 90^ confidence limits of the 
obsonred E value at tho times t^,^ t^, tg, and tg. The inference might be 
dravm that at these times energy loss processes, not considered in the fore- 
cast, are of importance. It is believed significant that a vrell defined 

^Additional computations of those data indicate that tho correct value of 
E is approximately idontioal with the forecast value. 
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aroa of oross-soa vras present at such location and tinos as to have inter- 
fered v/'ith the "svAsll” arriving at the times tg# t^, and tg from fetches A 
and B (Appendix III). Cross-seas are believed to cause the summation of the 
spectral components to become non-linear giving rise to the conversion of 
wave energy to tvirbulent onorgy by formation of whito-caps. Figure 18 
illustrates the location and orientation of the cross-sea relative to the 
generating fetches and the forecast sites at Davenport. Figure 19 displays 
the spectrum forecast for the cross-sea generating area at six hour inter- 

t 

vals over the period when it is most clearly defined. For reasons previously 
stated, computed spectra for the verification times under discussion, are 
not available. It is believed that such spectra would havo aided greatly in 
the dotonuination of the portion of the speotnim eiffooted by the cross-sea 
interference. However, it was noted that if the forecast swell energy vrere 
tobally erased in the frequency band of the cross-seas, energy reductions 
of the order of those observed in Figure 14 would have occurred. 

Figure 15 compares the forecast values of E with the observed values 
at site II, at a depth of 80 feet, for the times t^ through tg. 

By comparing the energy differences vdaich ocovir at the surface 
(Figure 14) with the energy difforoncos which occur at depth (Figure 15) 
it appears it might be permissible to draw inferences about the frequency 
ranges in which these differences occur. Those inferences aro based on 
the fact that the middle and lower frequency energy would bo less 
attenuated with depth by the hydro-dynamio filter than would the high 
frequency energy. Accordingly, if the relative error of the onorgy differ- 
onoos at the surface was of the same order as the relative error at depth, 
it might be inferred that the energy differences wore primarily in the 
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uiddlo and lair/ frequencies. By similar reasoning, if the relative error 
at depth v/as significantly less than the relative error at the surface, 
it might bo inferred that the energy differences vrere primarily in the high 
frequency end of the spectrum. If the relative error at depth is greater 
than the relative error at tho sea surface the possible causes are less 
obvious. A possible explanation might be that too much middle or low 
frequency energy and too little high frequency energy is forecast, the 
errors cancelling and giving close agroomont with the observed energy at 
tho surface. Then, as tho forecast spectrum is reduced by the hydro- 
dynamic filter, at depth tho surplus middle and low frequency energy which 
was forecast would bo relatively loss attenuated with depth than the high 
frequency energy. Then at depth, the forecast energy by virtue of the 
surplus middle frequency energy would shmv a greater relative error than 
at tho surface. 

It should be noted that it is quite possible that minor vrind systems 
such as sea-breeze circulations, may have developed between the six hour 
map intervals of this study which could have developed energy sources not 
taken into account. Further, because site II is some distance seavrard of 
site I, tho total energy at site II is greater as tho shoaling effect^ is 
less pronounced. Since, in the above discussion, the energy differences 
are compared bet^roon site I at the sea surface and site II at depth, the 
relative error as observed from Figure 15 will be someivhat greater than 
would be the relative error at depth at site I. 

Examining the relative errors of energy differences at the soa- 
surface and at depth for the times tj through tg, and applying the 



Refraction differences betv/eon site I and site II were assumed 
negligible. 



reasoning suggested above, the following comraents are offered; 

1* At the times t 3 and t 5 the relative erros at the sea-surfaoo and 
at depth are about the same. At the time t^ the relative error is much loss 
at depth than at the surface. Follav;ing tho reasoning suggested above it 
might be inferred that at times tj and tg important energy difforoncos 
oocurred in the middle frequencies and that at time t^ the energy difference 
T/as primarily in the high frequency range. As noted previously, it is 
believed that cross-sea interference acted on the "swell” arriving at these 
times. These comparisons give rise to the possible oonclusioi} that the 
forecast spectiatm ovorprediots in tho low and middle frequencies and imder- 
prodicts in the hi^ frequencies.' 

2. At the times tg and ty, the relative errer of the energy differ- 
ences was significantly greater at depth than at the surface. As suggested 
above, a possible explanation for this may bo an overforecast of energy in 
the middle and low frequencies and an underforeoast of energy in tho high 
frequency band. It may be noted -that these times are at or near an energy 
peak. It might bo suggested than that for higher wind speeds the assumed 
energy spectnm ovorprediots in the middle frequencies and underpredicts 

in the higher frequencies. At least such an assumption is consistent with 
those data, 

3, At tho times tg and tg the relative error at the surface is signi- 
ficantly larger algebraically than at depth. Following tho line of reason- 
ing used for time it might be inferred that important energy differences 
occurred in the high frequency band. The cause of this assumed higli 
frequency energy reduction is not readily apparent. It is noted, hovrover, 
that at both time t^ and tg tho observed energy is at a minimum, A 
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possible explanation might be that at lovror wind speeds the assumed energy 
spoctnun may ovorprediot in the high frequency end ef the spectrum* 

6D Summary 

In summary, it is the opinion of the authors that the forecast method 
for swell is eporationally useful for fetches of the typo active in this 
study* based on this investigation, it appears that the oiirves and pro- 
cedures of the method , vdion properly applied, vj-ill result in a vravo 
forecast of useful accuracy* Figure 22, a plot of the forecast significant 
wave height compared to the observed illustrates this conclusion* As noted 
in paragraph 6A above, the results are consistent with the assumption that 
energy reduction processes not considered in the forecast are at tinos 
important* In particular, it is believed that cross-sea intorferonce was 
responsible for significant onorgj'' reductions at certain verification times 
in this investigation. It is believed that a qualitative correction would 
be possible by reducing the forecast spoctrm by a fraction of the cross-sea 
spectrm developed* 

The data are also consistent with tho premise that, for increasing 
wind speeds, tho forecast energy spectrum is too high in the lovr and middle 
frequencies and possibly too Iovt in the high frequency band, and that for 
lofwor wind speeds, an opposite trend appears possible* 

6E Future Research 

IThon certain modifications in the computer programs discussed in 
Ghaptor 4 are completed, and accurate observational spectra are available, 
more dofinito conclusions will bo possible* It is believed that cross-sea 
effects may become quite clear from a comparison of tho forecast and tho 
observed spectra* It is tho opinion of the authors that extensions of tho 
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methods usod in this investigation can provide tho quality and quantity 
of data nooessary for tho evaluation of forecast methods and tho effects 
of attenuation* 
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FIGURE 2 WAVE SPECTRUM 
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TABLE Z (Continued) 
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-straight, C-cyclonic, A-anticyolonio 
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TABLE 3 (Continued) 
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APPENDIX I 



Tho purpose of this appendix is to tabulate values of [X (f)^^ in 
Equation (l). Chapter 1, as a function of v/ind spood and frequency* 



f (f)] 



(2rr)5f6 

Tjiiere C s 4*8 x 10^ cm^ sec”^ 



2g^ 

4Tr^f^2 



1-3 - 4*85 X 10~ 



The values of the ^(f)j^ for the solution of Equation (la) 



( 1 ) 

(la) 

are given in 



Table (8)* 
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APPENDIX II 



DESCRIPTION OF LOCALE AND FACILIirES AT 
THE FORECAST SITE, DAVENPORT, CALIFORNIA 

Tho University of California, vrith the pormission of Santa Cruz 
Portland Cement Company, planned, established, and maintained an ocean 
wave recording program at Davenport, California. As shown in Figure 15, 
the pier extended twenty-two hundred feet seaward of tho beach eind vms 
exposed to the large waves from tho North Pacific generating areas* 
Elaborate instrumentation and recording devices were used. Continuous 
ocean wave records from 8 Nov 1952 to 1 Dec 1954 vrere obtained from tho 
Beach Erosion Board (BEB), step-resistor gages (parallel type)£^15J and 
Mark IX CQ pressure heads. The locations of those t\To instruments are 
shown in Figure 13. The BEB was located at the seaward end of the pier* 
Depth immediately below the BEB was 50 feet* The Hark IX ims located in 
81^ feet of water (below MLn/). 

Fluctuations of the sea surface were transferred by electrical cir- 
cuitry to a recording center located on shore. Here the fluctuations were 
permanently recorded on Esterlino-Angus Continuous Record Charts.^ Two 
time scales wore used: 

6“ per 1 minute, called "fast time" 

6" per 1 hour, called "slov/' time" 

Fast time was recorded every twelve hours for tvrenty minutes beginning at 
0430Z and 1630Z of each day. tj_ (forecast time)was determined by the fast 
time records. 

A continual maintenance schedule was kept, which involved mainly 
cleaning and re-calibration. 

^See Figure 21 
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APPENDIX III 



DESCRIPTION OF THE SYNOPTIC SITUATION 
DURDIG AND PRECEDING THE FORECAST TBE 

As explained in Chapter 3, the forecast time began at 0430Z, 28 Jan 
1954 and ended at 0430Z, 1 Fob 1954. 

The synoptic situation affecting the forecast locality dviring and 
preceding the forecast time remained basically stationary. A low pressure 
center (average 1000 MB) existed over the Gulf of Alaska or British 
Colimbia from 1230Z, 22 Jan 1954 to 1230Z, 28 Jan 1954. A strong high 
pressure center (1040 MB) maintained itself over Alaska and the Boring Sea 
xTith a ridge extending southeastvrard down to 35 N latitude. Figure 7 
illustrates the general synoptic situation that existed before and remained 
nearly constant during the forecast period. 

The large cyclonic circulation in the Gulf of Alaska created a fetch 
of length (Flq.) 1320 nautical miles. Fetch width (P^) was 520 nautical 
miles \7ith the distance from the center of the leeward side to Davenport 
(R^) being 420 nautical miles. This fetch hereaPtor irill be referred to 
as Fetch ”A." Fetch A existed fully developed from 0030Z, 23 Jan 1954 to 
0630Z, 26 Jan 1954. 

At approximately 0630Z, 26 Jan 1954, due to a slight reorientation of 
the cyclonic circulation, F^q became 800 nautical miles with a correspond- 
ing increase in Rq to 960 naxitical miles. This change resulted in a 
redesignation to Fetch ”B.” All other geographic variables remained 
constant. Fetch B existed from 0630Z, 26 Jan 1954 to 1230Z, 28 Jan 1954. 

Data for fetches A and B v/ere recorded in Table 2, Chapter 3, and 
the fetches are illustrated in Figxire 7. 

Ylind is one of the most important variables in forecasting an energy 
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spoctrian. The general wind area remained fairly constant but wind varied 
before and during the forecast period in fetches A and D. The variation 
of the average wind speed with time in fetches A and B is illustrated in 
Figxire 17« 

Fig\ire 8 illustrates the general synoptic situation that existed from 
0030Z, 27 Jan 1954 to 1830Z, 27 Jan 1954. A fetch. Fetch C, was fitted 
to this circulation. Data resulting from this fetch are recorded in 
Table 2, Chapter 3. 

The wind field remained constant at 22 knots during the existence of 
Fetch C* 

Fetch A existed for 84 hours, B for 48 hours, and C for 18 hoTors. The 
subscript notation (A^,B^,Cj^) was'used to distinguish fetch sequence, with 
a map interval of 6 hours* The fetch life started when i = 1, and ended 
vdien i a 14, i » 8, and i » 4 respectively. 
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A study and evaluation of a 
method for forecasting ocean 
svrell waves using wave spectra. 




